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Abstract 

The influence of the HBT effect on the momentum spectra of in- 
dependently produced pions is studied using the method developed 
earlier for discussion of multiplicity distributions. It is shown that 
in this case all the spectra and multiparticle correlation functions are 
expressible in terms of one function of two momenta. It is also shown 
that at the critical point all pions are attracted into one quantum state 
and thus form a Bose-Einstein condensate. 



1. Introduction 

Several years ago, Pratt realized that the well-established phenomenon 
of HBT correlations @ can, under certain conditions, lead to Bose-Einstein 
condensation in multipion systems, which he called "pion laser". Since then 
the effect was investigated by several authors, including Pratt || |J (see || for 
an exhaustive list of references). Recently ||, employing the density matrix 
formalism, we discussed multiplicity distributions of independently emitted 
identical particles for arbitrary shapes of particle spectra in momentum and 
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in configuration space. The approach to the condensation point and the 
conditions for reaching it were investigated. In the present paper we extend 
the method of [[J to momentum spectra of independently produced identical 
pions. The explicit formula for the generating functional is written down and 
the resulting momentum spectra are discussed. 

Although the independent production mechanism is unlikely to be a re- 
alistic model of pion production, we feel that it is interesting to investigate 
its consequences. 

2. The basic ideas of the HBT effect, as applied to processes of particle 
production, were explained in 0, using the approach developed in [[/]]. The 
net result is the formula for the momentum distribution of n identical bosons 

%)==lEp (0) (?MP'). (1) 

where p^°\q,q') is the n-particle density matrix calculated ignoring the iden- 
tity of particle^ and normalized by the condition 

Tr\p®] = JdqpW(q,q) = l. (2) 

The sum extends over all permutations P and P' of particle momenta [qi, q n ] 

q- 

Consider now a system of n particles emitted independently. If we ignore 
the identity of particles, independent emission implies that the density matrix 
factorizes 

n 

P (0) (M') = IIP (0) fe.?D- (3) 

i=i 

Introducing this into ([!]) we have 

1 n 

n(«) = -EIIp (p) ((?p)i.(«p')*). ( 4 ) 

The information contained in @ is conveniently summarized in the form of 
the generating functional defined as 

11 ES*P'°>(n)W,.[l] 1 ' 



1 All quantities calculated with the identity of particles ignored will be called uncorrected 
and denoted by a superscript . 
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with 

W n [u] = J dqx...dq n Q(qi, q n )u(q 1 )...u(q n ). (6) 

Here p(°>(n) is the uncorrected multiplicity distribution and u(q) is an arbi- 
trary real nonnegative function of q. For u(q) = const = z the generating 
functional reduces to the generating function of the multiplicity distribution, 
§(z) discussed in [[|. 

The inclusive and exclusive distributions of n identical particles can be 
obtained by n-fold functional differentiation of the generating functional (|^) 
with respect to u(q) at u(q) = 1 and u(q) = 0, respectively Similarly 
n— fold differentiation of the logarithm of at u(q) = 1 gives the inclusive 
correlation function of the n-th order. 

To find the explicit expression for $[«] we observe that, given the formula 
(f|) for Q(q), for each permutation P of the momenta q±, ...,q n , the integral 
on the right hand side of @ factorizes into a product of contributions from 
all the cycles of P (as is well known, each permutation can be decomposed 
into cycles) . Let us denote the contribution from a cycle of length k by C k [u] . 
We have 

C k [u} = J d 3 q 1 ...d 3 q k u{q 1 )p {0) {q 1 ,q 2 )u{q 2 )p i0) {q 2 ,q 3 )...u{q k )p {0 \q k ,q 1 ) (7) 
The rest of the calculation is just combinatorics. 

We observe first that any two permutations which have identical par- 
titions into cycles give equal contributions. Let us consider the set of all 
permutations with a given partition into cycles. Denoting by n k the number 
of occurrences of a cycle of length k in the set of permutations considered, 
the contribution from all of them can be written as 

n< 1 » 

w<m = ntcM^p - »vr- x = «! n 4/-. (8) 

In the first equality the first factor is the integral, the second is the number 
of partitions of the n particles among the cycles, the third is the number of 
ways a cycle can be constructed from k particles and the last one corrects 
for the permutations of whole cycles. 

W n [w] is obtained by summing over partitions into cycles different 

from each other. 
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Until now we have considered a fixed multiplicity. As noted already in 
||, the sum over multiplicities can be explicitly performed if the uncorrected 
multiplicity distribution P™(n) is poissonian (as required for independent 
emission) 

. .71 

P®(n) = e- V %. (9) 
n! 

The result is an elegant formula for the generating functional fl5|): 

$M = exp I ^ ^ fc —J- ( 10 ) 

3. Let us now discuss the general properties of the particle distributions 
obtained from Eq.flIU|). 

The main result is that all inclusive distributions can be expressed in 
terms of a single function L(q, q'), defined as 

oo 

L(g,gO = E^V 0) ]W) ( n ) 

k=l 

where 

[p (0) f (<?,</) = I d 3 q 2 ...d 3 q kP ^(q,q 2 )p^(q 2 ,q 3 )...p^(q k ,q') (12) 



and the ^-dependence of L is not written explicitly. The single particle 
distribution is given by 

u(q)=L(q,q), (13) 

and the two-particle correlation function is 

K 2 (qi, q 2 ) = L(q x , q 2 )L(q 2 , q x ). (14) 

The general formula for the correlation functions reads 

K p (qi, ...q P ) = L(q u q 2 )L(q 2 , q 3 )....L(q p , q x ) 

+permutations of (q 2 , ,q p ). (15) 

It is not difficult to verify that by integrating fli~5|) over all momenta one 
recovers the formula for cumulants derived in || . 



These formulae represent a formidable constraint on the observed particle 
distributions: They basically say that all higher order correlations can be 
derived from the two particle correlation function. As they are valid for any 
model which assumes independent production, they were found by many 
authors in particular cases [§, [LC], [□]], for a full list, see ||. It is fair to say, 
however, that - since the independent production model is not expected to 
be a precise description of high energy interactions one expects violations 
of these relations, at least to some extent. Recent work by Eggers et al. 
[|T^] , points perhaps in this direction. The measurements of deviation of data 
from (|13|)-(|T5|) is of great interest, as it may indicate what is the dominant 
intraparticle correlation. 

Further discussion is greatly simplified if the matrix p^(q, q) is expressed 
in terms of its eigenvalues \ m and its eigenf unctions ip m (q). We haveQ 

P i0) (q,q') = J2M<lMM) (16) 

m 

and 

b (0) ]W) = E^(?)W t C(?')- (17) 

m 

Substituting this into (|ll]) and performing the sum over k we have 

L( qi ,q 2 ) = £VUfc)Va<fe)7^7r- (I 8 ) 

on 1 IS Am. 



It is clear from ([181) that L(g 1; g 2 ) and thus also all inclusive distributions 
become singular when i>X Q — ► 1 (Ao is the largest eigenvalue). In this limit, 
corresponding to Bose-Einstein condensation, L(q, q') is dominated by the 
first term in the sum and we have 

J ( \ ^(gi)^^^) , ff , ^ V 

L(qi,q2) = — r + L(q 1 ,q 2 ), (19) 

1 — UAo 

where L(qi,q2) remains bounded for u\q — > 1. Thus at the condensation 
point all the particles, except for a negligible fraction, are in the same state 
described by the eigenfunction ipo{q). Away from the condensation point, of 
course, the whole sum in (1181) must be carried out. 



2 We discuss here only the case of a discrete eigenvalue spectrum. The case of a 
continuous spectrum can be treated along the same lines H . 
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Let us now consider the important example, when particles are emitted 
in a pure state, i.e. 

p(q,q')=Hq)r(q')- (20) 
It follows from ([7[) that then 

C k [u] = (C 1 [u]) i (21) 
and thus the generating functional becomes 

One easily sees that this gives the momentum spectrum !?/'((/) | 2 independent 
of particle multiplicity and the geometric distribution of multiplicities 0. 

4. Let us now consider the uncorrected single particle density matrix 
of gaussian form, discussed already in several papers by Pratt [1], [3], [| and 
recently by other authors |5|, |TTJ. To simplify the notation, we shall re- 
strict ourselves to the one- dimensional problem. Generalization to the three- 
dimensional case is straightforward 0. We thus have 

where 

q + = ^(q + q'); q' = q-q'- (24) 

As easily seen, A 2 is the average value of the square of the particle momen- 
tum, and R 2 is the average value of the square of the space coordinate of the 
particle emission point. As is clear from the context, both A and R refer to 
the uncorrected distributions. The uncertainty principle implies that 

RA > \. (25) 

As explained in 0, the eigenf unctions of the density matrix ( p3|) are of 
the form 

ip m (q) = a m e-^ q2 H m (J —q), (26) 
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where H m (q) is the Hermite polynomial of order m and a m is the normalizing 



factor, given e.g. in [13| . 

The corresponding eigenvalues are 

A m = A (l - A ) m , m = 0,l,..., (27) 

where 

^Tnkm- 1 (28) 

is the greatest one. 

Eqs. (E5|) and (P7|) can be now used for the explicit calculation of L(q, q') 
[c.f. (|18D1 which, in turn, determines all particle distributions, as explained 
in the previous section. This is actually easier than it looks because the 
difficult sum over oscillating Hermite polynomials can be replaced by a sum 
over Gausssians with positive coefficients. To see this, we rewrite the formula 
^[181) as a double sum 



L(q, q') = £ 1, n (q)P m tf) £ i^mf • (29) 

m k=l 

Introducing (]27|) and reversing the order of summation we obtain 

L(q,q') = f: - ( ^ 0) " k PkM), (30) 
fe=l 1 — (1 — A J 



where 



2 _ (9 + ) 2 _1 p2 / --,2 



^ 5 ' )= l^J 6 ' (31) 

and Afc and i?fc are determined from the equations 



A fc A 21 -uj k \ 2RA + 1J 

The important lesson from this exercise is that even when the uncorrected 
distribution is described by a simple Gaussian, the resulting particle spectra 
are fairly complicated superpositions of an infinite number of Gaussians with 
varying width. 
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It follows from ([!2|) that Rk < R and < A for all k > 1. Consequently, 
the observed distributions are always narrower than the assumed uncorrected 
ones. 

Another interesting quantity is the average value < (q~) 2 > calculated 
from the two-particle correlation function, i.e., 

<(q-) 2 >=^-Jdq + K 2 (q, q ')(q-) 2 (33) 

where the cumulant K 2 is the integral of K 2 (q, q'). In the standard analysis 
of the data, 2 < (g~) 2 > is usually interpreted as an inverse of the average 
squared radius R e ff of the particle emission region. The ratio R 2 eff /R 2 cal- 
culated from (|33"D is plotted in Figure 1 for various values of RA. One sees 
that, at a fixed RA, R 2 ff decreases from R 2 to R 2 /2RA when < n > varies 
from to oo. Thus we conclude that, even fairly far from the critical point, 
the apparent size of the system R 2 tf (as determined from the two-particle 
correlation function) has little to do with the actual size of the system, given 
by R 2 . One sees also that, even at the fixed particle phase-space density 
< n > /RA, the effect substantially increases with increasing RA. 

5. Our conclusions can be summarized as follows. 

(i) It has been shown that the formulation of the problem of HBT inter- 
ference in terms of the density matrix provides an effective tool for discussion 
not only of particle multiplicities || but also of the momentum spectra. 

(ii) For independent particle production, the effects of HBT symmetriza- 
tion are expressed in terms of one function of two momenta. Consequently, 
all multiparticle correlations are expressible in terms of the single particle 
distribution and the two-particle correlation function. It would be very in- 
teresting to verify if this general relation holds - at least approximately - in 
the data. 

(iii) The symmetrized distributions reveal the existence of a critical point 
corresponding to the Bose-Einstein condensation, when almost all the par- 
ticles fall into one quantum state, corresponding to the lowest eigenvalue of 
the (uncorrected) density matrix. 

(iv) Due to symmetrization, the width of the momentum spectrum de- 
creases and the width of the two-particle correlation function increases. This 
effect becomes stronger, when the system approaches criticality. Close to the 
critical point the determination of the size of the emission region from the 
width of the correlation function is no longer possible. 
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(v) For a Gaussian uncorrected density matrix, the corrected momentum 
distribution is Gaussian both in the low density and in the high density limit, 
but non-Gaussian in the intermediate region. At low density the parameters 
A and R are determined by the widths of the momentum distribution and 
of the correlation function, respectively. At high density both these widths 
depend on the ratio 7c only. 
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Figure captions 

Figure 1. R e ff/R 2 plotted versus < n > for different values of RA. The 
marked points correspond to a fixed value of the particle phase-space density 
< n > /RA = 8. 
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